From hard to soft diffraction and return 
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The long standing mystery of smallness of diffractive dissociation of hadrons to large effective masses (the 
Pomeron-proton cross section is only 2mb) witnesses that the gluonic clouds of valence quarks are so small 
(ro = 0.3 fm) that soft interaction hardly resolves those gluons (diffraction is oc ro). A color-dipole light-cone 
(LC) approach is developed which incorporates a strong nonperturbative interaction of the LC gluons. The 
energy dependent part of the total hadronic cross section is calculated in a parameter-free way employing the 
nonperturbative LC wave functions of the quark-gluon Fock states. It rises with energy as s A , and we predict 
A = 0.17 ± 0.01, as well as the normalization. However, the energy independent part of the cross section related 
to inelastic collisions with no gluon radiated (gluons are not resolved) cannot be calculated reliably and we treat 
it as an adjustable parameter which is fixed fitting just one experimental point for total cross section. Then the 
energy dependence of the total cross section (the Pomeron part) and the elastic slope are fully predicted, as well 
as the effective Pomeron trajectory in impact parameter space, in a good agreement with data. These results 
naturally explain the ^-dependence of the diffractive DIS observed at HERA. Although diffraction is expected to 
be dominated by soft interactions the observed effective A is about twice as large as one (0.08) known for total 
cross sections. Diffractive excitations of large effective mass correspond to diffractive gluon radiation and should 
be associated with our calculated A. 



1. Diffraction: evidence for small gluonic 
spots in the proton 

The cross section of single diffraction in 
hadronic collisions (pp unless specified) can be ex- 
pressed in terms of the so called Pomeron-proton 
total cross section jj]J|, 



da(pp -> pX) Pp 



dt dM\ 



(1) 



What should one expect for af^f ? If to think 
about the Pomeron as a gluonic system its cross 
section contains the Cazimir factor 9/4 compared 
to the meson-nucleon cross section. Thus, one 
may expect er to | ) ^ 40 — 50 mb. Such a naive es- 
timate is, however, in a strict contradiction with 
data which show that the Pomeron interacts with 
a tiny cross section about 2 mb at high cm. en- 
ergy Mx This is the well know puzzling 
smallness of the triple-Pomeron coupling which 
was under intensive discussion back in 70s and is 
still waiting for a dynamical explanation. 



Although the unitarity corrections suppress 
diffraction they cannot bridge such a large gap 
between the expectation and data. Indeed, noth- 
ing surprising comes out from the analysis |J 
of the same data for the triple-Regge coupling 
IPIPR much larger that IP IP IP. The correspond- 
ing Reggeon-proton cross section af^ t « 20 mb 
agrees with the expectation, and the absorptive 
corrections do not disturb it much. 

The solution for the puzzling smallness of large- 
mass diffraction was suggested in Q . Due to color 
screening (color transparency) the Pomeron- 
proton cross section depends quadratically on the 
size of the gluonic system. Indeed, glueballs are 
expected to be bound tighter than quarkonia [^) , 
therefore to have smaller radiuses. Even the sim- 
ple relation || between the color triplet (octet) 
string tension and and the slope of Reggeon 
(Pomeron) trajectory, k = 1/(2-^0;^^), leads 
to the octet string tension which is at least four 
times larger than one for a triplet string (even 
more since the effective a' F is substantially in- 
creased by the unitarity corrections). 



2 



Thus, we should conclude that the strong non- 
perturbative interaction of gluons is responsible 
for smallness of large mass diffraction. 

It worth reminding that the Pomeron is not a 
particle, but a shadow of inelastic processes, and 
the concept of Pomeron-proton cross section is ill 
defined (except if the Pomeron is a true Regge 
pole). We use this concept here in order to im- 
pose a scale which shows that the diffractive cross 
section is abnormally small. In what follows we 
employ the light-cone (LC) QCD formalism which 
relates the observed suppression of diffraction to 
the smallness of gluonic spots inside the proton. 
Indeed, the observed dependence of the cross sec- 
tion Eq. (Q), 1/M|-, can only be due to diffractive 
radiation of gluons since they are vector particles. 
The interaction with the target should be suffi- 
ciently hard to resolve the gluonic fluctuation in 
a valence quark. The smaller is the gluon cloud 
of a valence quark, the more difficult is to ex- 
cite it. If the LC wave function of the proton 
is dominated by configurations in which gluons 
are located around the valence quarks at small 
transverse distances ~ tq, then the amplitude of 
diffraction is proportional to r§ and the cross sec- 
tion of diffractive gluon radiation is suppressed as 
Tq. Therefore, high-mass diffraction is a very sen- 
sitive tool to study the gluonic structure of the 
proton. The related formalism based on the LC 
Green function approach is presented below. 

2. From hard to soft diffraction 

The first correct pQCD calculation for large- 
mass diffraction cross section in DIS was per- 
formed in in the LC dipole representation. 
The diffractive cross section of gluon radiation 
q N — > q G N is given by Eqs. (64) of (| (see sim- 
ple derivation in Appendix A. 2 of H), 



M 2 da{qN -> qGN) 



81 



1024tt 



dM 2 dqi 



d 2 p 



q T =0 



* 'qG{a, P) (Tqq{p, s) 



,(2) 



C*<1 



where qx is the transverse momentum of the radi- 
ated gluon, M is the quark-gluon effective mass, 
a is the fraction of the initial LC momentum of 



the quark taken away by the gluon; dq q {p,s) is 
the cross section of interaction with a nucleon of 
a qq dipole with separation p and energy s. The 
perturbative quark-gluon (photon) LC wave func- 
tion is derived in ||] (see also in Using it one 
ends up with diffraction which grossly overesti- 
mates data for pp — > pX, as one might has ex- 
pected in view of discussion in previous section. 

Of course one should not apply pQCD methods 
to the soft reaction, at small qx the final quark 
and gluon cannot be treated as free particles since 
they should experience a strong nonperturbative 
interaction comoving with a large transverse sep- 
aration along the light cone. Propagation of such 
a qG pair from initial longitudinal coordinate z\ 
and separation p\ to final Z2 , P2 is described by 
the Schrodinger equation 0], 



% -7— G q G (zi , pi; Z2, p%) 
az2 



2pa (1 — a) 



+ Vg q (z 2 ,p 1 a) 



G 



qG(Zl, Pi] Z 2 , P2) 



(3) 



The first term in square brackets is the kinetic 
one, the LC potential is real for propagation in 
vacuum and describes the quark-gluon interac- 
tion. We choose it in the oscillator form in order 
to solve (||) analytically, 

6 4 (a)p 2 



ReV q G(z2,p,a) 



2pa(l — a) ' 



(4) 



where b 2 (a) — 6q + Abf a (1 — a). Since we are 
interested ina< 1, b(a) = b . 

As far as the Green function is known one can 
calculate the nonperturbative LC wave function 
for a quark-gluon Fock state, 



VqG&a) 



2tt pa(l — a) 



x / dz 1 [xTx)G qG (z 1 ,p 1 ;z2,p) 



Pi=0 



(5) 



where x, \ are the initial or final quark spinors, 
and the vertex function reads, 

r = i m q a 2 e* • (n x (?) + a e* ■ (er x V) 

- i(2-a)e*-V, (6) 
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which has a rather simple form, 

2i la s t ■ r 

y 



* gG (f,a < 1) 



(7) 



where e is the polarization vector of the gluon. 
This wave function recovers the perturbative one 
H in the limit of bo — > 0. 

Comparing the cross section of diffractive gluon 
radiation with data one arrives at a rather large 
value of b — 0.65 GeV corresponding to a 
short mean quark-gluon separation tq — 1/bg = 
0.3 fm. This is one of the central results of 

3. Small gluonic spots and the total pp 
cross section 

Using the same approach one can calculate the 
cross section of gluon bremsstrahlung (nondiffrac- 
tive) by a qq meson 0,O,pT| , 



hN 
1 

9 
4 



= / dOlr, 



d 2 R 



o 

da . 
— / d~r 

a 



Vq G (R + r,a) 



x ag (R + r) + * ?G (r , a) ag (r) 
- Re^* qG (r,a)^g G (R + r,a) 



agJR + f)+agJr)-ag Q (R) 



(8) 



where R and r are the qq and Gq transverse 
separations respectively. Making use of relation 
(r 2 ) <C (R 2 ) one can keep only the second term 
in the curly brackets what makes further calcula- 
tions rather simple. 

The cross section of gluon bremsstrahlung by a 
quark interacting with a nucleon is found pOjhTl 
to be, 



1 



rad 



n=l 





' s N 




n 9 
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37T ' 


v s o, 
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Crl 



(9) 



Since the mean separation ro is small, one can use 
the approximation Oqq(r) w Cr 2 , where C ~ 2.3 
can be evaluated perturbatively. The QCD cou- 
pling a s was averaged over the radiation spec- 
trum with a result a s — 0.38 — 0.43 in a good 



agreement with the critical value a c = 0.43 eval- 
uated in Uj. 

To obtain the total pp cross section one should 
add the cross section of pp interaction without 
any gluon radiation. This contribution, (Tq P is 
pure nonperturbative and cannot be evaluated re- 
liably. We treat it as a free parameter. Thus, the 
total pp cross section takes the form, 

A 



with 



27 



Cr 2 



4 a 

A = — - = 0.17 ±0.01 , 

6TT 



(10) 



(11) 



, -<-9Cr 2 /4. 
This value of A is about twice as large as the 
one suggested by the data for the energy depen- 
dence of total pp cross sections when the simple 
parameterization af^ t oc s A is applied. However, 
the radiative part is a rather small fraction of the 
total cross section (at medium high energies) and 
the presence of the energy independent term <7o 
substantially reduces the effective A e //. 

4. Unitarization and comparison with data 

The cross section (|Io| ) apparently violates the 
Froissart bound and one should introduce unitar- 
ity corrections ||l3|| . Unfortunately, this is not a 
well defined procedure since different recipes can 
be found in the literature. 

The simplest known way to restore unitarity is 
to eikonalize the partial amplitude "fp(b, s), 



Im r p (b, s) = 1 — exp — Im (6, s) 



(12) 



At very high s this amplitude approaches the 
black disk limit ^,lmT P (s,b) -► S[R 2 (s)-b 2 ] , 
with radius, R(s) = roAln(s/so). Correspond- 
ingly, at asymptotic energies, 



Aln 



> 



(13) 



all hadronic cross sections reach the maximal 
universal energy growth allowed by Froissart- 
Martin's bound. 



hN 



(s) -> 2nr 2 A 2 In 2 



(14) 
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The eikonalization procedure (12) would be 
suitable if the the incoming hadrons were eigen- 
states of the interaction JL4J. Hadrons, how- 
ever, are subject to diffractive off-diagonal ex- 
citation, and the eikonal form of unitarization 
should be corrected in a way similar to Gribov's 
inelastic corrections for hadron-nucleus cross sec- 
tions. The lowest order unitarity correction in 
dl2] ) comes from the quadratic term in the ex- 
ponent expansion of T(b, s). It has to be modi- 
fied using the AGK cutting rules to include single 
diffraction, 



Imlp = 



1 



= Im 72P - - ( Im 7ip 



1 + D(s) 



(15) 



where D(s) ~ cr s d(s)/cr e i(s) is approximately 0.25 
in the ISR energy range and decreases slightly 
with energy oc s -0 04 |15| due to stronger unitar- 
ity corrections. Asymptotically, as s — * oc, D(s) 
vanishes since cr e ;(s) oc ln 2 s and <r s d{s) oc Ins. 

The inelastic corrections to higher order terms 
in the expansion ( |l2] ) are poorly known. A simple 
way to keep (15) and to include diffraction into 
the higher terms is to modify (02) as, 



ImT]p(b, s 
|l — exp 



1 



- 1 



f £>(*) 

£>(«)) Im 7JP (M)]} (16) 

which is known as quasi-eikonal model [pj. 

Following Jl6| we assume that the ^dependence 
of the pp elastic amplitude is given by the Dirac 
electromagnetic formfactor squared. Correspond- 
ingly, the mean square radius (r 2 h ) evaluated in 
Jig is smaller than (r 2 h ). 

For the dipole parameterization of the formfac- 
tor the partial elastic amplitude which is related 
via unitarity to the nth term in the radiation cross 
section (ph takes the form, 



(17) 



where K^{y) is the third order modified Bessel 
function and y = b^8/B n . The slope param- 
eter grows linearly with n due to random walk 
of radiated gluons with a step 7q in the impact 
parameter plane, B n = 2(f 2 h }/3 + nr 2 . 



In order to calculate the total cross section, 
otot = 2/ d 2 blmT(b, s), one needs to fix the 
energy independent term with n = in (|l0|). 
This can be done comparing with the data for 
atot at any energy sufficiently high to neglect 
Reggeons. We use the most precise data 17 at 
= 546 GeV and fix a = 39.7 mb. 

The predicted energy dependence of a^ t shown 
by the dashed curve in Fig. 1 is in good agreement 
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Figure 1. Data for total pp (full circles) and pp 
(open circles) cross section and the prediction 
of Eq. ( ]l7| ) for the energy dependence of the 
Pomeron part (dashed curve) and including 
Reggeon contribution (solid curves. 



with the data at high energies, but apparently 
needs Reggeon corrections towards low energies. 

In order to improve description of data at 
medium high energies (ISR) we added the con- 
tribution of leading Reggeons parameterized in 
the standard way with common parameters for 
pp and pp, the intercept fixed at or(Q) = 1/2 
and the slope parameter a' R = 0.9 GeV~ 2 . The 
Reggeons are added directly in the partial elastic 
amplitude, 

Iml(s,b) = ImT]p(s,b) 

+ r R (s,b) \l - lmT F (s,b)} , (18) 
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with a proper screening by absorptive corrections. 
The results are shown by the solid curves of Fig. 1 
(pp bottom curve and pp upper curve). 

Employing Eq. ( |l6| ) we can also predict the 
slope of elastic scattering at t = 0, B e i(s) = 
(b 2 ) /2, where averaging is weighted by the partial 
amplitude (|l^). The results exhibit good agree- 
ment when compared with the pp and pp data in 
Fig. 2. 
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Figure 2. Data for the elastic slope and our 
predictions. The upper and bottom curves 
and open and full circles correspond to pp and 
pp respectively. 



Note that the the slope essentially depends on 
our choice of (f 2 h ) in (|l7|), however the predicted 
energy dependence, i.e. the effective value 
is fully defined by the parameter ro fixed in 
Indeed, each radiated gluon makes a small "step" 
~ r 2 — 0.1 fm 2 in the impact parameter plane 
leading to the rising energy dependence of the 
elastic slope. Eventually, at very high energies the 
approximation of small gluon clouds should break 
down. Nevertheless, the mean number of gluons 
in a quark (n) = A ln(s/so) remains quite small in 
the energy range of colliders. It is only (n) = 0.7— 
1 at the ISR and reaches about two gluons at the 
Tevatron. Correspondingly, the mean square of 



the quark radius grows from 0.06 fm 2 to 0.18 fm 2 
which is still rather small compared to the mean 
square of the charge radius of the proton. 

5. Pomeron trajectory in impact parame- 
ter space 

Actually, we know more about the elastic par- 
tial amplitude than just integral characteristics 
like the total cross section and slope, we know 
the shape of 6-dependence. To compare the b- 
distribution and its development with energy di- 
rectly with experimental data an analysis similar 
to one which has been done by Amaldi and Schu- 
bert jl8| was performed recently in [^TJ . The data 
for differential elastic cross section was fitted, in- 
cluding both the real and imaginary parts of the 
amplitude, and then Fourier transformed to im- 
pact parameter representation. As different from 
Pq| , no model dependence was involved (geomet- 
rical scaling was assumed in fllg|| ) and data from 
the SppS at yfs = 540 GeV were included in ad- 
dition to data from the ISR. The data for partial 
elastic amplitude as function of impact parameter 
are plotted in Fig. 3 in comparison with our pre- 
dictions based on Eq. (|l8|). Agreement between 
the data and our predictions is remarkably good, 
especially if to recall that the Pomeron part has 
no free parameters, except one, <r , adjusted to 
the total cross section measured at one energy 
y/s = 546 GeV 0. Both the predicted shape of 
the partial amplitude and its energy development 
are confirmed by the data. 

Addition of the SppS data at y/s = 540 GeV 
into the analysis plays a crucial role, it substan- 
tially increases the energy range important for 
study of energy dependent effects. In particu- 
lar it helps to extract from the data the effective 
Pomeron trajectory ap(b) = 1 + A(b) as func- 
tion of impact parameter, describing energy de- 
pendence s A ( b ) of the partial elastic amplitude. 
The results are depicted in Fig. 4. The solid 
curves show the slope of energy dependence cal- 
culated with Eq. Jl8| ) without any adjustment. 
The agreement is rather good. 

The simple parameterization of the amplitude 
by a single Pomeron pole without unitarity cor- 
rections fitted to data for total cross section and 
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Figure 3. The imaginary part of the par- 
tial amplitude Imr(6) as function of impact 
parameter at different cm. energies. The 
solid curves are the theoretical prediction with 
Eq. and the dashed curves show the con- 
tribution of Reggeons. 
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Figure 4. The effective Pomeron trajectory 
in impact parameter space. Black and open 
points correspond to different parameteriza- 
tion of the t-dependent elastic amplitude. Our 
prediction and of the single Pomeron pole 
model are shown by solid and dashed curves 
respectively. 



elastic slope |19(] corresponds to the dashed curve. 
Although this parameterization is very successful 
describing total cross section it poorly reproduces 
energy dependence of the partial amplitude at dif- 
ferent impact parameters. At the same time we 
should notice that the wide spread prejudice that 
this model leads to a precocious breaking of uni- 
tarity at small b is incorrect. Wee see from Fig. 4 
that A (6 = 0) predicted by this simple model is 
rather small. In the case of geometrical scaling it 
would be zero. 

6. Back to hard diffraction: Do we under- 
stand the results from HERA? 

Diffraction of highly virtual photons, 7* p — > 
Xp observed at HERA is expected to be dom- 
inated by soft processes pQ] . Indeed, the domi- 
nant hard fluctuations of the photon, have a small 



transverse size ~ l/Q 2 , therefore the diffraction 
cross section is tiny, ~ 1/Q 4 - On the other 
hand, soft fluctuations of large size which appear 
very rarely with probability <~ 1/Q 2 (for trans- 
verse photons), interact with a large hadronic 
cross section, therefore they dominate the "hard" 
DIS diffraction. Consequently, one might expect 
that the cross section should rise with energy (at 
fixed Q 2 ) like the hadronic cross sections, oc s 08 . 
However, a much steeper slope of energy depen- 
dence was observed at HERA. According to the 
analysis in ^lj data from the HI and ZEUS ex- 
periments lead to 

A = 0.2 ±0.02 (HI) , 

A = 0.13 ±0.04 {ZEUS) . (19) 

The source of excitations with large effective 
masses is diffractive gluon bremsstrahlung related 
via unitarity to the energy dependent second term 
in the elastic amplitude (total cross section) in 
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Eq. ([L0|) . At the same time, small mass diffraction 
is governed by excitation of the qq fluctuation in 
the photon without gluon radiation (yet it may in- 
clude radiation of energetic gluons) . This diffrac- 
tive component is related by unitarity to the first 
energy independent term in Eq. (|To|). There- 
fore, for large excitation masses we should expect 
about the same value of A given by Eq. (pT[), 
rather than the averaged effective value ~ 0.08 
fitted to the total cross section. 

Formulating differently, we can say that the 
softest part of interaction responsible for the con- 
stant term in the total cross section Eq. ( [To| ) is 
related to the large size R of the valence quark 
skeleton. Such interaction is too soft to resolve 
the gluon clouds of small size Tq. It needs a 
semi-hard interaction with transverse momentum 
~ 1 GeV to shake off the gluons from the quarks. 
This part of interaction is responsible for the en- 
ergy dependent term in (|l0|) and for diffractive 
gluon radiation, i.e. excitation of large effective 
mass. 

7. Summary 

Our main observation is existence of small glu- 
onic clouds surrounding the valence quarks in 
light hadrons. The most soft part of interaction 
with transverse wave length longer than these 
spots cannot either resolve them or shake them 
off the quarks. As an example, it can be treated 
in the string model as a crossing and rearrange- 
ment of the hadronic strings. The corresponding 
part of the total cross section is energy indepen- 
dent since the size of the valence quark skeleton 
does not grow with energy. 

On the contrary, the gluon spots do rise with 
energy since the gluon population increases and 
they perform a Braunian motion in the transverse 
plane. The gluon bremsstrahlung cross section 
also rises since each gluon can be radiated at dif- 
ferent rapidities with about the same probability. 
This is the source of energy dependence of the 
total cross section. 

Such a two-scale picture of interaction leads to 
presence of the two terms in the total cross section 
Eq. (fl0|). The first one is energy independent, and 
the second one rises as power of energy with quite 



a large exponent given by (p. l|). This second term 
is fully predicted, not only its energy dependence, 
but also the normalization. Only the first energy- 
independent term is treated as a free parameter 
and is fixed by fitting just one experimental point. 

Our calculations reproduce well the energy de- 
pendence of the total pp cross section, and the 
elastic slope. Moreover, the data for differential 
elastic scattering Fourier transformed to impact 
parameter representation also well agree with our 
predictions, including the effective Pomeron tra- 
jectory in impact parameter space. 

The two components of the total cross sec- 
tion can be well separated in diffractive dissoci- 
ation. Indeed, the small mass diffraction corre- 
sponds to excitation of the quark ensemble with- 
out gluon radiation, i.e. it is energy independent. 
Only diffractive excitation of gluon clouds leads 
to the large mass dissociation, called in Regge 
phenomenology triplc-Pomcron diffraction. This 
part is rather small, but steeply rises with energy 
with the power given by (|ic|). This indeed has 
been observed by experiments at HERA. 

We have found a new regime where one can 
perform calculations employing the smallness of 
the size tq = 0.3 fm of the gluonic spots com- 
pared to large interquark separation R ~ 1 fm. 
On the contrary, in DIS one deals with an oppo- 
site limit when the gluon clouds are much larger 
than the interquark separation. In that case the 
cross section can be estimated using pcrturbativc 
QCD provided that Q 2 > Ql (Q ~ 1 GeV), and 
s > so (so ~ IGeV 2 ), but x = Q 2 / s < 1. De- 
pendent on the approximation used, two models 
for the Pomeron are known, the BFKL [^2) and 
double-leading-log approximation. The total vir- 
tual photoabsorption cross section is predicted to 
rise steeply with energy as is confirmed by data 
from HERA. 

Our conclusion about smallness of gluonic 
spots in light hadrons is in a good accord with 
what has been already observed in models for 
nonperturbative QCD vacuum. It has been found 
employing the QCD sum rules j23| that the glu- 
onic formfactor of the proton has very weak Q 2 
dependence corresponding to a small radius of the 
gluon distribution ~ 0.3 fm. A small gluon corre- 
lation radius ^0.3 fm also appears in the lattice 
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calculations [|24| . It is also predicted by the model 
of instanton liquid |2E}| , Experimental observa- 
tion of a small cross section for large mass soft 
diffractive dissociation which has led the analyses 
Q to a small value of rg can be treated more gen- 
erally, in particular as a confirmation for a small 
size cloud of any kind of gluonic vacuum fluctu- 
ations surrounding the valence quarks. They are 
usually referred to as constituent quarks. 

Note that the two-scale structure of the energy- 
dependent total cross section has to exist also in 
other models attempting to incorporate the non- 
perturbative effects [^6 27 1. 
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